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Gravitation as a fundamental interaction that governs all phenomena at large and very small scales, but still not 
well understood at a quantum level, is a cardinal missing link in unification of all physical interactions. Discovery 
of the present acceleration of the Universe, the dark matter and dark energy problems are also a great challenge 
to modern physics, which may lead to a new revolution in it. Integrable multidimensional models of gravitation 
and cosmology make up one of the proper approaches to studying basic issues and strong field objects, the early 
and present Universe and black hole (BH) physics in particular [TJE1E]. Our main results within this approach are 
described for both cosmology and BH physics. Problems of the absolute G measurements and its possible time and 
range variations are reflections of the unification problem. 

The choice, nature, classification and precision of determination of the fundamental physical constants as well 
as their role in a transition, expected in 2011, to new definitions of the main SI units, supposed to be based on 
fundamental physical constants and stable quantum phenomena, are described. The problem of temporal variations 
of constants is also discussed, temporal and spatial variations of G in particular. A need for further absolute 
measurements of G, its possible range and time variations is pointed out. The multipurpose space project SEE is 
briefly described, aimed at measuring G and its stability in space and time, with precision 3-4 orders better than at 
present. It may answer many important questions posed by gravitation, cosmology and unified theories. A project 
of a laboratory experiment for testing possible deviations from the Newton law is also presented. 
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1. Introduction 

Studies in the previous century in the field of gravi- 
tation were mainly devoted to theoretical studies and 
experimental verification of general relativity (GR) and 
alternative theories of gravity with a strong stress laid 
on relations between macro- and micro-world phenom- 
ena or, in other words, between classical gravitation and 
quantum physics. Very intensive investigations in these 
fields were done in Russia by M.A. Markov, K.P. Sta- 
niukovich, Ya.B. Zeldovich, A.D. Sakharov and their 
colleagues starting from mid-60s. As a motivation there 
were: the existence of singularities in cosmology and 
black hole physics, the role of gravity at large and very 
small (Planckian) scales, attempts to create a quantum 
theory of gravity as for other physical fields, problem 
of possible variations of fundamental physical constants 
etc. A lot of work has been done in such areas as [3]: 

— exact solutions with different fields as sources in GR, 

— particlclikc solutions with a gravitational field, 

— quantum field theory in a classical gravitational 

background, 

— quantum cosmology with fields (e.g., scalar), with 

the cosmological constant etc., 

— self-consistent treatment of quantum effects in cos- 

mology, 

— development of alternative theories of gravity: scalar- 

tensor, gauge, with torsion, bimetric etc., 
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— possible variations of fundamental physical constants 

As our main results of this period, one can men- 
tion [4] the first quantum cosmological model with a cos- 
mological constant (creation from nothing) (1972); the 
first classical cosmological models for a conformal scalar 
field (1968) and quantum cosmological models with min- 
imal and conformal scalar fields (1971), the first nonsin- 
gular cosmological model with spontaneous symmetry 
breaking of a nonlinear conformal scalar field (1978-79), 
exact solutions for nonlinear electrodynamics, including 
the Born-Infcld one, the first exact solution for dilaton- 
type interaction with an electromagnetic field in GR, the 
first nonsingular particlelike purely field solution with 
gravity (1979). One should also mention the conclusion 
that only G may vary with respect to atomic system 
of measurements, or in the Jordan-Brans-Dicke frame 
(1978). 

Since all attempts to quantize general relativity in 
a usual manner failed and it was proved that it is non- 
renormalizable, it became clear that a promising trend 
goes along the lines of unification of all physical interac- 
tions which has started in the 70s. About this time, ex- 
perimental studies of gravity in strong fields and gravita- 
tional waves started, giving a powerful speed-up in the- 
oretical studies of such objects as pulsars, black holes, 
QSO's, active galactic nuclei, the early Universe etc., 
which continue now. 

But nowadays, when we think of the most important 
lines of future developments in physics, we can foresee 
that gravity will be essential not only by itself, but as 
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a missing cardinal link of some theory unifying all ex- 
isting physical interactions: gravitational, weak, strong 
and electromagnetic ones. Even in experimental ac- 
tivities, some crucial next generation gravitational ex- 
periments verifying predictions of unified schemes will 
be important. Among them are: MICROSCOPE and 
STEP — testing the cornerstone Equivalence Princi- 
ple, SEE — testing the inverse square law (or new non- 
Newtonian interactions), testing possible time variations 
of the Newtonian constant G , measurement of the ab- 
solute value of G with unprecedented accuracy [TOUTl"] . 
All these experiments became tests of not only Grav- 
ity itself, but unified models of physical interactions as 
well. Of course, the gravitational-wave problem, veri- 
fication of torsional, rotational (GPB), 2nd order and 
strong field effects remain important as well. 

We can also predict that the studies of gravity it- 
self and within the unified models will give, in the next 
century and millennium, even more applications to our 
everyday life, as the electromagnetic theory gave us in 
the 20th century after very abstract fundamental inves- 
tigations of Faraday, Maxwell, Poincare, Einstein and 
others, who had never dreamed of such enormous appli- 
cations of their works. 

Another very important feature which may be envis- 
aged is an increasing role of fundamental physics studies, 
gravitation, cosmology and astrophysics in particular, in 
space experiments [12 • Unique micro-gravity environ- 
ments and the modern technology outbreak give nearly 
a perfect place for gravitational experiments which suf- 
fer a lot on Earth from its relatively strong gravitational 
field and gravitational fields of nearby objects due to the 
fact that there is no way of screening gravity. 

In the development of relativistic gravitation and 
dynamical cosmology after A. Einstein and A. Fried- 
mann, we may notice three distinct stages: first, studies 
of models with matter sources in the form of a perfect 
fluid, as was originally done by Einstein and Friedmann. 
Second, studies of models with sources as different phys- 
ical fields, starting from the electromagnetic and scalar 
ones, in both classical and quantum setting (see [4]). 
And third, which is really topical now, application of 
ideas and results of unified models for treating funda- 
mental problems of cosmology and black hole physics, 
especially in high energy regimes and for explanation of 
the greatest challenge to modern physics, explaining the 
present acceleration of the Universe, the so-called dark 
energy problem. Multidimensional gravitational models 
play an essential role in the latter approach. 

The necessity of studying multidimensional models 
of gravitation and cosmology [TJ[2] is motivated by sev- 
eral reasons. First, the main trend of modern physics 
is unification of all known fundamental physical interac- 
tions: electromagnetic, weak, strong and gravitational 
ones. In the recent decades, there has been a significant 
progress in unifying weak and electromagnetic interac- 
tions and some more modest achievements in GUT, su- 



persymmetric, string and superstring theories. 

Now, theories with membranes, p-branes and more 
vague M-theory are being created and studied. Having 
no definite successful theory of unification now, it is de- 
sirable to study the common features of these theories 
and their applications to solving basic problems of mod- 
ern gravity and cosmology. Moreover, if we really be- 
lieve in unified theories, the early stages of the Universe 
evolution and black hole physics, as unique superhigh 
energy regions, and possibly even the low energy stage, 
when we observe the present acceleration, are the most 
proper and natural arena for them. 

Second, multidimensional gravitational models, as 
well as scalar-tensor theories of gravity, are theoretical 
frameworks for describing possible temporal and range 
variations of fundamental physical constants [4, fjJElIZ]- 
These ideas have originated from the earlier papers by 
E. Milne (1935) and P. Dirac (1937) on relationships be- 
tween the phenomena of micro- and macro- worlds, and 
up till now they are under thorough study both theoret- 
ically and experimentally. The possible discovery of the 
fine structure constant variations is now under critical 
further investigation. 

Lastly, applying multidimensional gravitational mod- 
els to basic problems of modern cosmology and black 
hole physics, we hope to find answers to such long- 
standing problems as singular or nonsingular initial 
states, creation of the Universe, creation of matter 
and its entropy, the cosmological constant and coin- 
cidence problem, origin of inflation and specific scalar 
fields which may be necessary for its realization, the 
isotropization and graceful exit problems, stability and 
nature of fundamental constants [5lfT2l[T3]. the possible 
number of extra dimensions, their stable compactifica- 
tion, new revolutionary data on present acceleration of 
the Universe, dark matter and dark energy etc. 

Bearing in mind that multidimensional gravitational 
models are certain generalizations of GR which is tested 
reliably for weak fields up to 0.0001 and partly in strong 
fields (binary pulsars), it is quite natural to inquire 
about their possible observational or experimental win- 
dows. From what we already know, among these win- 
dows are: 

— possible deviations from the Newton and Coulomb 

laws, or new interactions, 

— possible variations of the effective gravitational con- 

stant with a time rate smaller than the Hubble 
one, 

— possible existence of monopole modes in gravita- 

tional waves, 

— different behaviour of strong field objects, such as 

multidimensional black holes, wormholes and p- 
branes, 

— standard cosmological tests, 

— possible nonconservation of energy in strong field 

objects and accelerators if the brane-world ideas 
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about gravity in the bulk turn out to be true etc. 

Since modern cosmology has already become a unique 
laboratory for testing the standard unified models of 
physical interactions at energies far beyond the level 
of existing and future manmade accelerators and other 
installations on Earth, there exists a possibility of using 
cosmological and astrophysical data for discriminating 
between future unified schemes. Data on possible time 
variations or deviations from the Newton law as a new 
important test should also contribute to the unified 
theory choice. 

As there exist no accepted unified models, in our ap- 
proach [Tlf2 |l 1 411X5] we have adopted simple (but general 
from the viewpoint of the number of dimensions) mod- 
els, based on multidimensional Einstein equations with 
or without sources of different nature: 

— the cosmological constant, 

— perfect and viscous fluids, 

— scalar and electromagnetic fields, 

— their possible interactions, 

— dilaton and moduli fields with or without potentials, 

— fields of antisymmetric forms (related to p-branes) 

etc. 

The main objective of our programme was and is to ob- 
tain exact self-consistent solutions (integrable models) 
for these systems and then to analyze them in cosmo- 
logical, spherically and axially symmetric cases. In our 
view, this is a natural and most reliable way of studying 
highly nonlinear systems. It is done mainly within Ric- 
mannian geometry. Some simple models in integrable 
Weyl geometry and with torsion were studied as well. In 
many cases, we tried to single out models which do not 
contradict the available experimental or observational 
data on variations of G . In some cases we have used the 
methods worked out for arbitrary dimensions in study- 
ing 4D models as well. 

As our model [2(5] , we use n Einstein spaces of con- 
stant curvature with sources as (m + 1) -component per- 
fect fluid (or form fields, ), cosmological or spherically 
symmetric metrics, in manifolds obtained as a direct 
product of factor spaces of arbitrary dimensions. Then, 
in the harmonic time guage, we show that the Einstein 
multidimensional equations are equivalent to Lagrange 
equations with a nondiagonal (in general) minisuper- 
space metric and some exponential potential. After pos- 
sible diagonalization of this metric, we perform reduc- 
tion to sigma-model and Toda-like systems and further 
to Liouville, Abel, generalized Emden-Fowler equations 
etc. and try to find exact solutions. We suppose that the 
behaviour of extra spaces is the following: they can be 
constant or dynamically compactified (e.g., toroidally), 
or large, but with barriers, walls etc. 

So, we have been realizing the programme in arbi- 
trary dimensions (from 1988) [UHI3KMI2] • 



In cosmology: 

we have obtained exact general solutions of multidimen- 
sional Einstein equations with sources: 

— A, A + scalar field (e.g. nonsingular, dynamically 

compactified, inflationary); 

— perfect fluid (PF), PF + scalar field (e.g. nonsingu- 

lar, inflationary solutions); 

— viscous fluid (e.g. nonsingular, generation of mass 

and entropy, quintessence and coincidence in a 2- 
component model); 

— stochastic behaviour near the singularity, billiards in 

Lobachevsky space, that D = 11 is critical and <p 
destroys billiards (1994). 

For all the above cases, Ricci-flat solutions were ob- 
tained for any n , and, in addition, solutions with cur- 
vature in one factor space; with curvatures in 2 factor 
spaces, solutions are known only for total D = 10, 11; 

— fields: scalars, dilatons, forms of arbitrary rank 

(1998) — inflationary solutions, A generation by 
forms (p-branes) [71] : 

— the first billiards for dilaton-forms (p-branes) inter- 

action (1999); 

— quantum systems (solutions of the WDW equation 

[IB] ) for all the above cases where classical solu- 
tions were obtained; 

— dilatonic fields with potentials, billiard behaviour for 

them as well. 

For many of these integrable models, we have calculated 
the time variation of the effective gravitational constant 
and compared it with the present experimental bounds, 
which allowed choosing particular models or singling out 
some classes of solutions. 

Solutions depending on r in any dimensions: 

— generalized Schwarzschild and Tanghcrlini (BHs are 

singled out), solutions with a minimally coupled 
scalar field ip (no BHs); 

— generalized Reissner-Nordstrom (BHs are also sin- 

gled out), the same plus ip (no BHs); 

— multitcmporal solutions; 

— for dilaton-like interaction of ip and electromagnetic 

fields (BHs exist only in a special case); 

— stability studies for the above solutions (only BH 

ones are stable); 

— the same for dilaton-forms interaction (p-branes); 

stability was found only in some cases, e.g., for 
a single form. 

The PPN parameters for most of the models were cal- 
culated. 

Theory of experiments: 

Space and laboratory experiments aimed at testing 
a possible violation of the Newton law and raising the 
precision of the absolute value of the Newton constant 
G determination were suggested and worked out. 
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2. Multidimensional Models 

The history of the multidimensional approach begins 
with the well-known papers of T.K. Kaluza and O. Klein 
on 5-dimensional theories, which aroused an interest in 
studies of multidimensional gravity. These ideas were 
continued by P. Jordan who suggested to consider the 
more general case c/55 ^ const leading to a theory with 
an additional scalar field. They were, in some sense, a 
source of inspiration for C. Brans and R.H. Dicke in their 
well-known work on a scalar-tensor gravitational theory. 
After their work, a lot of investigations have been per- 
formed with material or fundamental scalar fields, both 
conformal and non-conformal (see details in [3]). 

A revival of ideas of many dimensions started in 
the 70s and is continuing now, entirely due to the de- 
velopment of unified theories. In the 70s, an inter- 
est in multidimensional gravitational models was mainly 
stimulated by (i) the ideas of gauge theories leading to 
a non-Abelian generalization of the Kaluza-Klein ap- 
proach and (ii) by supergravitational theories. In the 
80s, the supergravitational theories were "replaced" by 
superstring models. Now, it is heated by expectations 
connected with the overall M-theory. In all these theo- 
ries, 4-dimensional gravitational models with extra fields 
were obtained from some multidimensional model by 
dimensional reduction based on decomposition of the 
manifold as 

M = M 4 x M int , (1) 

where M 4 is a 4-dimensional manifold and M;„t is some 
internal manifold (mostly considered to be compact). 

The earlier papers on multidimensional gravity and 
cosmology dealt with multidimensional Einstein equa- 
tions and with a block-diagonal cosmological or spher- 
ically symmetric metric defined on the manifold M = 
Rx M x ... x M„ of the form 

n 

g = -dt<Z>dt + ^2a 2 r (t)g r , (2) 

r=Q 

where (M r ,g r ) are Einstein spaces, g r is a metric on 
M r , r = 0, ... ,n. In some of them, simple scalar fields 
and a cosmological constant were also included [16] , 

Such models are usually reduced to pseudo- Euclidean 
Toda-like systems with the Lagrangian 

L=-G ij x i V-'£A k e<* i (3) 

k=l 

and the zero-energy constraint E = . 

It should be noted that pseudo-Euclidean Toda-like 
systems are not yet well studied. There exists a special 
class of equations of state that gives rise to Euclidean 
Toda models [17]. 

Cosmological solutions are closely related to solu- 
tions with spherical symmetry [IS]. The first multidi- 
mensional generalization of such type was considered 



by D. Kramer and rediscovered by A.I. Legkii, D.J. 
Gross and M.J. Perry (and also by Davidson and Owen). 
Moreover, the scheme of obtaining the latter is very sim- 
ilar to the cosmological approach [HUH]. We continued 
these investigations in detail. In [5D], the Schwarzschild 
solution was generalized to the case of n internal Ricci- 
flat spaces, and it was shown that a black hole configura- 
tion takes place when the scale factors of internal spaces 
are constants. It was also shown there that a minimally 
coupled scalar field is incompatible with the existence of 
black holes. In [3T], an analogous generalization of the 
Tanghcrlini solution was obtained, and an investigation 
of singularities was performed in [25] . These solutions 
were also generalized to the electrovacuum case with and 
without a scalar field [22,23. 24;. Here, it was also proved 
that BHs exist only when the scalar field is switched off. 
Deviations from the Newton and Coulomb laws were 
obtained, depending on the mass, charge and number 
of dimensions. In [24 , spherically symmetric solutions 
were obtained for a system of scalar and electromagnetic 
fields with a dilaton-type interaction, and also devia- 
tions from the Coulomb law were calculated, depending 
on the charge, mass, number of dimensions and the dila- 
tonic coupling. Multidimensional dilatonic black holes 
were singled out. A theorem was proved in [24j that 
"cuts" all non-black-hole configurations as being unsta- 
ble under even monopole perturbations. In [26], the ex- 
tremely charged dilatonic black hole solution was gener- 
alized to multi-centre (Majumdar-Papapetrou) configu- 
rations when the cosmological constant is non-zero. 

We note that for D = 4 the pioneering Majumdar- 
Papapetrou-type solutions with a conformal scalar field 
and an electromagnetic field were considered in [27j . 

At present, there exists a special interest in the 
so-called M-theory etc. These theories are "super- 
membrane" analogues of the superstring models in 
D = 11 and higher dimensions. The low-energy limit 
of these theories leads to multidimensional models with 
p-branes. 

Exact solutions with "branes" 

In our papers, several classes of the exact solutions were 
considered for the multidimensional gravitational model 
governed by the Lagrangian 

C = R[g] -2A- h aP g MN d M ¥ a d N ^ 

Tia • 

a 

Here g is a metric, F a = dA a are forms of ranks n a , 
tp a are scalar fields, and A is a cosmological constant 
(the matrix h a p is invertible). 

Supergravities. For certain field contents with dis- 
tinguished values of the total dimension D , ranks n a , 
dilatonic couplings A Q and A = , such Lagrangians ap- 
pear as "truncated" bosonic sectors (i.e. without Chern- 
Simons terms) of certain supergravitational theories or 
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a low-energy limit of superstring models [128,29. For 
D = 11 supergravity [25] (considered now as a low- 
energy limit of a conjectured M -theory |30]). we have a 
metric and a 4-form in the bosonic sector. For D = 10, 
one may consider type I supergravity with a metric, a 
scalar field and a 3-form; type IIA supergravity, with 
bosonic fields of type I supergravity, called the Neveu- 
Schwarz-Neveu-Schwarz (NS-NS) sector, and addition- 
ally a 2-form and a 4-form Ramond-Ramond (R-R) sec- 
tor; type IIB supergravity with bosonic fields of type 
I supergravity (the NS-NS sector) and additionally a 1- 
form, a 3-form and a (self-dual) 5-form (the R-R sector). 
It is now believed that all five string theories (I, IIA, IIB 
and two heterotic ones with gauge groups G — Eg x E& 
and Spin(32)/Z2) [29] as well as 11-dimensional super- 
gravity [28] are limiting cases of M-theory. All these 
theories are conjectured to be related by a set of dual- 
ity transformations: S— , T— (and more general U— ) 
dualities. 

It was proposed earlier that the IIB string may have 
originate in a 12D theory, known as F-theory (Vafa). 
A low energy effective (bosonic) Lagrangian for F- 
theory was also suggested. The field content of this 
12-dimensional field model is the following: a metric, 
a scalar field (with a negative kinetic term), a 4-form 
and a 5-form. In our work |31| . a chain of so-called 
-Bp-models in dimensions D = 11,12,... was sug- 
gested. The Brj -model contains I = D — 11 scalar 
fields with negative kinetic terms (i.e. the so-called 
"phantom" fields) coupled to (I + 1) different forms of 
ranks 4, . . . , 4 + 1 . These models were constructed using 
p-brane intersection rules that will be discussed below. 
For D = 11 (I — 0), the Br> -model coincides with 
the truncated bosonic sector of D = 11 supergravity. 
For D — 12 (I = 1) , it coincides with the truncated 
D = 12 model. We have conjectured in [31] that these 
Br> -models for D > 12 may correspond to low-energy 
limits of some unknown Fr> -theories (analogues of M— 
and F -theories). 

Description of the models. In our review |T5], cer- 
tain classes of p-brane solutions to field equations corre- 
sponding to the Lagrangian (4), obtained by us earlier, 
were presented. 

These solutions have block-diagonal metrics defined 
on a D -dimensional product manifold, i.e., 

n 

9 = c 2 V + e2</ V > M o x Mi x ... x M„, (5) 

i=l 

where g° is a metric on Mq (our space) and g l are fixed 
Ricci-flat (or Einstein) metrics on Mi (internal spaces, 
i > 0). The moduli 7, (jf and the scalar fields ip a are 
functions on Mq , and fields of forms are also governed 
by several scalar functions on Mq . Any F a is supposed 
to be a sum of monoms, corresponding to electric or 
magnetic p-branes (p-dimensional analogues of mem- 
branes), i.e., the so-called composite p-brane ansatz is 



considered [3H[33] . (In the non-composite case we have 
no more than one monomial for each F a . ) p = corre- 
sponds to a particle, p = 1 to a string, p — 2 to a mem- 
brane etc. The p-brane world volume (world line for p = 
0, world surface for p — 1 etc.) is isomorphic to some 
product of internal manifolds: Mr = M ix x . . . x M ik 
where 1 < i\ < . . . < < n and has the dimension 
p + 1 = d n + . . . + d lk = d(I) , where J = {i u . . . , i k } is 
a multi- index describing the location of the p-brane, and 
di = dim Mi. Any p-brane is described by the triplet 
(p-brane index) s = (a, v, I) , where a is the colour in- 
dex labelling the form F a , v = e(lectric) , m(agnetic) 
and / is the multi-index defined above. For the elec- 
tric and magnetic branes corresponding to a form F a , 
the world volume dimensions are d(I) = n a — 1 and 
d(I) — D — n a — 1 , respectively. A sum of these dimen- 
sions is D — 2. For D — 11 supergravity we get d(I) = 3 
and d(I) = 6, corresponding to an electric M2-brane 
and a magnetic M5-brane. 

Sigma-model representation. In our paper [34] . 
the model under consideration was reduced to a grav- 
itating self-interacting sigma-modcl with certain con- 
straints. The sigma-model representation for the non- 
composite electric case was obtained earlier in [321133] . 
for the electric composite case see also |35j). 

The cr-model Lagrangian, obtained from (4), has the 
form [34] 

£ CT = R[g a ] - G AB g^d^ A d„o- B 

-J2?s exp(-2C/ s ).g°^^<i> s ^$ s - 2V, (6) 

s 

where (a A ) = (0\</? Q ), V is a potential, (Gab) are 
components of the (truncated) target space metric, e s = 
±1, 

U s = U s A a A = ditf ~ XsX asa ^ a 

are linear functions, $ s are scalar functions on Mo (cor- 
responding to the forms), and s — (a s ,v s , I s ) . The pa- 
rameter Xs = +1 for the electric brane (v s = e) and 
Xs — — 1 for the magnetic one (v s = m). 

The pure gravitational sector of the sigma model was 
considered earlier in our paper [18] . For p-brane appli- 
cations, g° is Euclidean, (Gab) is positive-definite (for 
do > 2 ) , and e s = — 1 if pseudo-Euclidean (electric and 
magnetic) p-branes in a pseudo-Euclidean space-time 
are considered. The sigma model (6) may also be con- 
sidered for the pseudo-Euclidean metric g° of signature 
(—,+,...,+) (e.g. in studies of gravitational waves). 
In this case, for a positive-definite matrix (Gab) and 
e s = 1 , we get non-negative kinetic energy terms. 

The brane U -vectors. The co- vectors II s play a 
key role in studying the integrability of the field equa- 
tions [331 SD] and possible existence of stochastic be- 
haviour near the singularity, see our paper 50 . An 
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important mathematical characteristic here is the ma- 
trix of scalar products (U S ,U S ) — G AB U A U% , where 
(G AB ) — {Gab)^ 1 ■ The scalar products for the co- 
vectors U 8 were calculated in [M] (for the electric case 
see [32E135]): 

(uw) = d(i s ni s ,) + d ^M 

where {h al) ) = (7i aj g) -1 ; s = (a s , v s , I s ) , s' = 
(a s <, ty , I s >) . They depend on brane intersections (the 
first term), the brane world- volume dimensions and 
the total dimension D (the second term), the scalar 
products of dilatonic coupling vectors and the electro- 
magnetic types of branes (the third term). As will be 
shown below, the so-called "intersections rules" (i.e., re- 
lations for d(I s (~1 I s ')) are determined by the scalar 
products of II s -vectors. 

Cosmological and spherically symmetric 
solutions 

A family of general cosmological-type p-brane solutions 
with n Ricci-flat internal spaces was considered in our 
paper [42], where a generalization to the case of n — 
1 Ricci-fiat spaces and one Einstein space of non-zero 
curvature (say, Mi ) was also obtained. These solutions 
are defined up to solutions to Toda-type equations and 
may be obtained using the Lagrange dynamics following 
from our sigma-model approach 3 lj . The solutions from 
[42 contain a subclass of spherically symmetric solutions 
(for Mi = S dl ) . Special solutions with orthogonal and 
block-orthogonal [35] sets of U -vectors were considered 
previously in our works [31] and [T4l|42], respectively. 
(For the non-composite case, see [37 ] l38 ] ) and references 
therein.) 

Toda solutions. In [3T], a reduction of p-brane cosmologio; 
type solutions to Toda-like systems was first performed. 
General classes of p-brane (cosmological and spherically 
symmetric) solutions related to Euclidean Toda lattices 
associated with Lie algebras (mainly A m , C m ones) 

were obtained by us in [3HES1I1S1I1710H] ' 

A class of space-like brane (S -brane) solutions (re- 
lated to Toda-type systems) with a product of Ricci-flat 
internal spaces and S'-brane solutions with special or- 
thogonal intersection rules were considered in [581159] . 
and solutions with accelerated expansion (e.g. with a 
power-law and exponential behaviour of scale factors) 
were singled out. 

Black-brane solutions. In our papers [47]I48|. a fam- 
ily of spherically symmetric solutions was investigated, 
and a subclass of black-hole configurations related to 
Toda-type equations under certain asymptotical condi- 
tions was singled out. These black-hole solutions are 
governed by the functions H s (z) > defined on the 
interval (0,(2/i) _1 ), where /i > is the extremality 



parameter, and obey the set of differential equations 
(equivalent to Toda-type ones) 

with the following boundary conditions: 

(i) i? s ((2/i)- 1 -0) = H sQ E (0, +oo); (ii) H s (+0) = 1, 

s G S. Here B s ^ and (A ss >) is a quasi-Cartan 
matrix. 

In Refs. [46,47,48] , the following hypothesis was put 
forward: the functions H s are polynomials when the 
intersection rules 

d(l s n i s .) = _ XsXs ,\ asa \ asl0 h^ 

+ ±(U S ' 1 U S ')A SS ,, s^s'. 

correspond to semi-simple Lie algebras, i.e., when (A ss ') 
is a Cartan matrix. Here, 

( 2(u;u') \ 

s, s' G S , is a quasi-Cartan matrix. 

This hypothesis was verified for Lie algebras: A m , 
C m +i , m = 1,2, . .., in |47 , 48J . It was also confirmed 
by special black-hole "block orthogonal" solutions con- 
sidered earlier in [T4"Il4l5] . 

In our papers [15] 07] [35], explicit formulae for the 
solution corresponding to the algebra A2 are presented. 
These formulae are illustrated by two examples of A2 - 
dyon solutions: a dyon in D = 11 supergravity (with 
M2- and M5-branes intersecting at a point) and a 
Kaluza-Klein dyon. Extremal configurations (e.g., with 
a multi-black-hole extension) were also obtained. 
I- We note that special black-hole solutions with or- 
thogonal [/-vectors were considered in [36j (for the non- 
composite case) and [31]. These solutions have ana- 
logues in models with a multi-component perfect fluid 
[UE155]. 

The black-brane solution, corresponding to the Lie 
algebras C2 and A3 were obtained in [52"] . 

In [36], some propositions related to i) interconnec- 
tion between the Hawking temperature and the singu- 
larity behaviour, and ii) multi-temporal configurations 
were proved. 

It should be noted that the polynomial structure was 
also found for the so-called flux-brane solutions which 
occur as generalizations of the well-known Melvin solu- 
tion. 

Cosmological models in diverse dimensions. Scalar 
fields play an essential role in modern cosmology though 
the problem of their origin still exists. They are at- 
tributed to inflation models of the early Universe and to 
models describing the present stage of the accelerated 
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expansion. There is no unique candidate potential of the 
minimally coupled scalar field. Typically a potential is 
taken as a sum of exponentials. Such potentials appear 
quite generically in a large class of theories: Kaluza- 
Klcin models, supergravity and string/M-theories. 

A single exponential potential was extensively stud- 
ied in the Friedmann-Robertson- Walker (FRW) 4D 
model containing both a minimally coupled scalar field 
and a perfect fluid with a linear barotropic equation of 
state. The attention was mainly focussed on the quali- 
tative behaviour of solutions, stability of exceptional so- 
lutions to curvature and shear perturbations and their 
possible applications within the known cosmological 
scenarios such as inflation and scaling ("tracking"). In 
particular, it was found by phase plane analysis that for 
"flat" positive potentials there exists a unique late-time 
attractor in the form of a scalar-dominated solution. It 
is stable within homogeneous and isotropic models with 
non-zero spatial curvature with respect to spatial cur- 
vature perturbations and provides power-law inflation. 
For "intermediate" positive potentials, a unique late- 
time attractor is the scaling solution where the scalar 
field "mimics" a perfect fluid, adopting its equation of 
state. The energy density of the scalar field scales with 
that of the perfect fluid. Using our methods for mul- 
tidimensional cosmology, the problem of integrability 
by quadratures of the model in 4 dimensions was also 
studied. Four classes of general solutions, when the 
parameter characterizing the steepness of the potential 
and the barotropic parameter obey some relations, were 
found [66 . For the case of multiple exponential poten- 
tial of the scalar field plus dust, an integrable model in 
4D was obtained in [67] . 

As to scalar fields with a multiple exponential po- 
tential in any dimensions, it is not yet well studied, al- 
though a wide class of exact solutions were obtained in 
our papers [53] 69]. In our recent work [57], the be- 
haviour of this system near the singularity was studied 
using a billiard approach suggested earlier in our pa- 
pers [56U50] . A number of S-brane solutions were found 
in [HHH!]. See details for 2-component D-dimensional 
integrable models in Refs. [70 l l64 ] l65] ) . 

Quite a different model with a dilaton, branes and a 
cosmological constant, with static internal spaces, was 
investigated in [71], where possible generation of the 
effective cosmological constant by branes was demon- 
strated. A model with variable equations of state was 
found in [68] . with acceleration in our space and com- 
pactification of internal spaces. 

Cosmological models with time variations of G. 

As was mentioned before, cosmological models in scalar- 
tensor and multidimensional theories are frameworks for 
describing possible time variations of fundamental phys- 
ical constants due to scalar fields present explicitly in 
STT and present initially or/and generated by extra 
dimensions in multidimensional approach. In [74], we 
have obtained solutions for a system of conformal scalar 



and gravitational fields in 4D and calculated the present 
possible relative variation of G at the level of less than 
10~ 12 /year. Later, in the framework of a multidimen- 
sional model with a perfect fluid and 2 factor spaces (our 
3D space could be open, closed or flat) and an internal 
6D Ricci-flat space, we obtained the same limit for such 
variation of G [S]. 

We have also estimated the possible variations of the 
gravitational constant G in the framework of a general- 
ized (Bergmann-Wagoner-Nordtvedt) scalar-tensor the- 
ory of gravity on the basis of the field equations, without 
using their special solutions. Specific estimates were es- 
sentially related to the values of other cosmological pa- 
rameters (the Hubble and acceleration parameters, the 
dark matter density etc.), but the values of G/G com- 
patible with modern observations 10 -12 /year [90] were 
not exceeded. 

In [55], we continued the studies of models in arbi- 
trary dimensions and obtained relations for G in a mul- 
tidimensional model with Ricci-flat internal space and 
a multi-component perfect fluid. A two-component ex- 
ample (dust + 5-brane) was considered. It was shown 
that G/G is less than 10~ 12 /year. Expressions for G 
were also considered in a multidimensional model with 
an Einstein internal space and a multicomponent perfect 
fluid [92] • In the case of two factor spaces with non-zero 
curvatures without matter, a mechanism for prediction 
of small G was suggested. The result was compared 
with exact (l+3+6)-dimensional solutions which we ob- 
tained earlier. 

A multidimensional cosmological model describing 
the dynamics of n + 1 Ricci-flat factor-spaces in 
the presence of a one-component anisotropic fluid was 
considered in [ 116] . The pressures in all spaces were sup- 
posed to be proportional to the density: pi — Wip, i — 
0, ...,71. Solutions with an accelerated power-law expan- 
sion of our 3-space Mo and small enough variation of 
the gravitational constant G were found. These solu- 
tions exist for two branches of the parameter w$ . The 
first branch describes superstiff matter with wo > 1 , the 
second one may contain phantom matter with wq < — 1 , 
e.g., when G grows with time, so this branch can de- 
scribe not the present epoch but rather earlier stages. 

Similar exact solutions, but nonsingular and with an 
exponential behaviour of the scale factors, were consid- 
ered in |117j for the same multidimensional cosmological 
model describing the dynamics of n + 1 Ricci-flat factor 
spaces Mi in the presence of a one-component perfect 
fluid. Solutions with accelerated exponential expansion 
of our 3-space Mq and small variation of G were also 
found. 

Exact S-brane solutions with two electric branes and 
two phantom scalar fields on the manifold 

M = R+ x K x Mi x M 3 x M 4 x M 5 . (7) 

were obtained and studied in [118] . We obtained asymp- 
totic accelerated expansion of our 3D factor space and 
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variations of G obeying the present experimental con- 
straints, G/G < 10- 12 /year. 

A D-dimensional cosmological model with several 
scalar fields and an antisymmetric (p + 2) -form was 
also considered |119j . For dimensions D — Am + 1 = 
5, 9, 13, ... and p — 2m — 1 = 1, 3, 5, ... , we obtained a 
family of new cosmological type solutions with a 4m- 
dimensional oriented Ricci-flat submanifold N of Eu- 
clidean signature. These solutions are characterized by 
a self-dual or anti-self-dual parallel charge density form 
Q of rank 2m defined on N. The (sub)manifold N may 
be chosen to be either Kahler or hyper-Kahler, or an 8- 
dimensional manifold of Spin(7) holonomy. A general- 
ization of these solutions to a chain of extra (marginal) 
Ricci-flat factor spaces was also presented. Solutions 
with accelerated expansion of extra factor spaces were 
singled out. Certain examples of new solutions for IIA 
supergravity and for a chain of Bp -models in dimen- 
sions D = 14, 15, ... were considered. 

Spherically symmetric solutions, black holes and 
PPN parameters. In [34], it was shown that, af- 
ter dimensional reduction on the manifold Mq x Mi x 
... x M n and when the composite p-brane ansatz is 
considered, the problem is reduced to the gravitating 
self-interacting a model with certain constraints. For 
electric p-branes, see also [33J [35] ( m [33, the com- 
posite electric case was considered). This representa- 
tion may be considered as a powerful tool for obtain- 
ing different solutions with intersecting p-branes. In 
[3H2D], Majumdar-Papapetrou type solutions were ob- 
tained (for the non-composite electric case see [33J and 
for the composite electric case see [35]). These solutions 
correspond to Ricci-flat (Mi,g l ) , i = 1, . . . , n and were 
generalized to the case of Einstein internal spaces [31] . 
These solutions take place when certain orthogonality re- 
lations (on the couplings parameters, brane dimensions 
and the total dimension) are imposed. In this situation, 
a class of cosmological and spherically symmetric solu- 
tions was obtained 31j. Solutions with a horizon were 
considered in detail in [3M3T] . 

It should be noted that multidimensional and multi- 
temporal generalizations of the Schwarzschild and Tan- 
ghcrlini solutions were considered in [53J[75], where the 
generalized Newton formulae in the multitemporal case 
were obtained. 

We have also calculated the Post-Newtonian Param- 
eters (3 and 7 (the Eddington parameters) for gen- 
eral spherically symmetric solutions and black holes in 
particular [14 . These parameters, depending on p- 
brane charges, their world-volume dimensions, dilaton 
couplings and the total number of dimensions may be 
useful for physical applications. 

Some specific models in classical and quantum mul- 
tidimensional cases with p-branes were analyzed in [311 
[37] , Exact solutions for a system of scalar fields and 
fields of forms with dilatonic type interactions for gen- 
eralized intersection rules were studied in [47] , where the 



PPN parameters were also calculated. Other problems 
connected with observations were studied in [601163] , and 
general properties of BHs and wormholes in a brane 
world in [62MI]. 

Stability analysis for solutions with p-branes was 
carried out in [35J85J. It was shown there that, for 
some simple p-brane systems, multidimensional black 
branes are stable under monopolc perturbations while 
other (non-BH) spherically symmetric solutions turned 
out to be unstable. 

Below we mainly dwell upon some problems of fun- 
damental physical constants, the gravitational constant 
in particular, upon the SEE and laboratory projects for 
measuring G and its possible variations and briefly on 
some theoretical models with variations of the effective 
gravitational constant. 

3. Fundamental physical constants 

1. In any physical theory we meet constants which 
characterize the stability properties of different types of 
matter: objects, processes, classes of processes and so 
on. Some of them cannot be calculated via other phys- 
ical constants. These constants are important because 
they arise independently in different situations and have 
the same value, at any rate within accuracies we have 
gained nowadays. That is why such constants are called 
the fundamental physical constants (FPCs) [4, 12J . It is 
impossible to define strictly this notion. It is because 
the constants, mainly dimensional ones, are present in 
certain physical theories. In the process of scientific 
progress, some theories are replaced by more general 
ones with their own constants, and there arise relations 
between old and new constants. So, we cannot speak of 
an absolute choice of FPC, but rather only of a choice 
corresponding to the present state of the physical sci- 
ences. 

Really, before creation of the electro-weak interac- 
tion theory and some Grand Unification Models, the 
following choice of FPCs was considered: 

c, H, a, G F , g si m p (or m e ), G, H, p, A, k, /, (8) 

where a, Gf, g s and G are constants of the electro- 
magnetic, weak, strong and gravitational interactions, 
H , p and A are cosmological parameters (the Hubble 
constant, the mean density of the Universe and the os- 
mological constant), k and / are the Boltzmann con- 
stant and the mechanical equivalent of heat which play 
the role of conversion factors between temperature on 
the one hand, energy and mechanical units on the other. 
After adoption (in 1983) of a new definition of the me- 
ter (A = ct or I = ct), this role is also partly played 
by the speed of light c. It is now also a conversion 
factor between units of time (frequency) and length, it 
is defined with absolute accuracy, i.e., zero uncertainty 
(with the new suggested definitions of basic units of the 
International System of Units (SI), such a role may also 
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be played by ft and Na , where Na is the Avogadro 
number [51] . 

Now, when the theory of electro-weak interactions 
has a firm experimental basis and we have some good 
models of strong interactions, a more preferable choice 
is as follows: 

ft, (c), e, m e , 6 W , Gf, 6 c , Aqcd, 

G, H, p, A, k, I (9) 

and maybe the three Kobayashi-Maskawa angles 62 , ^3 
and 6. Here 9 W is the Weinberg angle, 6 C is the Cabbibo 
angle and Aqcd is the cut-off parameter of quantum 
chromodynamics. Of course, if a theory of the four in- 
teractions known now will be created (M-theory or some 
other), then we will probably have another choice. As wc 
see, the macroscopic constants remain the same, though 
in some unified models, i.e. in multidimensional ones, 
they may be related in some manner (see below). From 
the point of view of these unified models, the above- 
mentioned ones are low-energy constants. 

All these constants are known with different accu- 
racies. The most precisely defined constant was and 
remains to be the speed of light c : its uncertainty was 
10~ 10 while now it is defined with an absolute accuracy. 
The atomic constants e, ft, m and others are deter- 
mined with errors 10 -6 ~ 10~ 8 , G up to 10~ 4 or even 
worse, 6 W up to 10~ 3 ; the accuracy of H is about a 
few per cent. Other cosmological parameters (FPCs): 
mean density estimations vary within 2 per cent; for A 
we now have data that its corresponding energy den- 
sity exceeds the matter density (0.7 vs. 0.3 of the total 
Universe mass). Here are some recent estimates from 
observational cosmology: 

— Total mean density : 0.98 < fi to t < 1-08. 

— Today's Hubble parameter: H = 0.72 ± 0.07. 

— Dark energy density parameter: Qde = 0.7. 

— For dark matter: SIdm — 0.26. 

— For baryonic matter: fig = 0.04. 

— For radiation: CIr = 5x 10 -5 . 

— Power spectrum index: n = 0.970 ± 0.023. 

— Equation of state coefficient: w = pj p < —0.78. 

As to the nature of the FPCs, we can mention several 
approaches. One of the first hypotheses belongs to J. A. 
Wheeler: in each cycle of the Universe evolution, the 
FPCs arise anew along with physical laws which govern 
this evolution. Thus the nature of the FPC and physical 
laws are connected with the origin and evolution of our 
Universe. 

A less global approach to the nature of the dimen- 
sional constants suggests that they are needed to make 
physical relations dimensionless or they are measures of 
asymptotic states. Really, the speed of light appears in 
relativistic theories in factors like v/c, at the same time 
velocities of usual bodies are smaller than c , so it also 



plays the role of an asymptotic limit. The same sense 
have some other FPCs: ft is the minimal quantum of 
action, e is the minimal observable charge (if we do not 
take into account quarks which are not observable in 
free state) etc. 

Finally, FPCs or their combinations may be consid- 
ered as natural scales determining the basic units. While 
earlier the basic units were chosen more or less arbitrar- 
ily, i.e., the second, metre and kilogram, now the first 
two are based on stable (quantum) phenomena. Their 
stability is believed to be ensured by the physical laws 
which include FPCs. There appeared similar sugges- 
tions for a new reproducible realization of a kg, to fix 
the values of Na or other constants, e.g. ft [99] . 

Another interesting problem, which is under discus- 
sion, is why the FPCs have values in a very narrow range 
necessary for supporting life (stability of atoms, stellar 
lifetime etc.). There exist several possible explanations. 
First, that it is a good luck, no matter how improbable 
is the set of FPCs. Second, that life may exist in other 
forms and with another FPC set, of which we do not 
know. Third, that all possibilities for FPC sets exist in 
some universe. And the last but not least: that there is 
some cosmic fine tuning of FPCs, some unknown physi- 
cal processes bringing them to their present values in a 
long-time evolution, cycles etc. 

An exact knowledge of FPCs and precision measure- 
ments are necessary for testing the main physical theo- 
ries, extension of our knowledge of nature and, in a long 
run, for practical applications of fundamental theories. 
Within this, certain theoretical problems arise: 

1) development of models for confrontation of theory 

with experiment in critical situations (i.e., for ver- 
ification of GR, QED, QCD, GUT or other unified 
models); 

2) setting limits on spacial and temporal variations of 

the FPCs. It is becoming especially important 
now, with the idea to introduce new basic units 
of the International System of Units (SI), based 
completely on fundamental physical constants. 

Of course, raising the precision of their absolute values 
is always a permanent task. 

As to a classification of the FPCs, we can set them 
into four groups according to their generality: 

1) Universal constants, such as ft which divides all 

phenomena into quantum and non-quantum ones 
(micro- and macro- worlds) and to a certain extent 
c, which divides all motions into relativistic and 
non-relativistic ones; 

2) coupling constants like a, 6 W , Aqcd and G; 

3) constants of elementary constituencies of matter like 

m e , m w , m x , etc., and 

4) transformation multipliers such as k , / and partly c 

(conversion from the second to the metre). Soon 
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there may be more after modernization of SI in 
2011 — the values of h and Na may be fixed 
with zero uncertainty. 

Of course, this division into classes is not absolute. 
Many constants move from one class to another. For 
example, e was a charge of a particular object, the elec- 
tron, class 3, then it became a characteristic of class 
2 (electromagnetic interaction, a — e 2 /(he) in combi- 
nation with h and c); the speed of light c has been 
in nearly all classes: from 3 it moved into 1, then also 
into 4. Some of the constants ceased to be fundamental 
(i.e. densities, magnetic moments, etc.) as they are 
calculated via other FPCs. 

As to the number of FPCs, there are two opposite 
tendencies: the number of "old" FPCs is usually dimin- 
ishing when a new, more general theory is created, but 
at the same time new fields of science arise, new pro- 
cesses are discovered in which new constants appear. So, 
in the long run, we may come to some minimal choice 
which is characterized by one or several FPCs, maybe 
connected with the so-called Planck parameters — com- 
binations of c, h and G (the natural, or Planck system 
of units [HE]): 

L = (hG/c 3 ) 1/2 ~ 10~ 33 cm, 
m L = (ch/2G) 1/2 ~ 1(T 5 g, 

t l = L/c ~ 10~ 43 s. (10) 

The role of these parameters is important since mi 
characterizes the energy of unification of the four known 
fundamental interactions: strong, weak, electromag- 
netic and gravitational ones, and L is a scale where 
the classical notions of space-time lose their meaning. 
There are other ideas about the final number of FPC 
(2, 1, or none, L. Okun' et al.) Of course, everything 
will depend on a future unified theory. 

2. The problem of the gravitational constant G mea- 
surement and its stability is part of a rapidly devel- 
oping field, called gravitational-relativistic metrology 
(CRM). It has appeared due to the growth of measure- 
ment technology precision, spread of measurements over 
large scales and a tendency to unification of fundamen- 
tal physical interactions, where the main problems arise 
and are concentrated on the gravitational interaction. 
This was first formulated in [7]. 
The main subjects of GRM are: 

— general-relativistic models on different astronomical 

scales: Earth, the Solar system, galaxies, cluster 
of galaxies, cosmology; 

— time transfer, VLBI, space dynamics, relativistic as- 

trometry etc. (pioneering works were done in Rus- 
sia by Arifov and Kadyev and Brumberg in the 
60s); 

— development of generalized gravitational theories 

and unified models for testing their effects in ex- 
periments; 



— measurement of fundamental physical constants, G 

in particular, and their stability in space and time; 
MICROSCOPE, STEP, SEE... 

— fundamental cosmological parameters as fundamen- 

tal constants: cosmological models studies (quint- 
essence, k-essence, phantom, multidimensional 
models), measurements and observations; projects 
PLANCK, ... 

— gravitational waves (detectors, sources...); LIGO, 

VIRGO, TAMA, LISA, RADIOASTRON,... 

— basic standards (clocks) and other modern preci- 

sion devices (atomic and neutron interferometry, 
atomic force spectroscopy etc.) in fundamental 
gravitational experiments, especially in space for 
testing GR and other theories: rotational, tor- 
sional and second-order effects (need an uncer- 
tainty of 10~ 6 or better), e.g. LAGEOS, Gravity 
Probe B, ASTROD, LATOR etc. 

We are now at the level of 2.3 • 10~ 5 in measuring the 
PPN-parameter 7 and 5 • 10~ 4 for j3 ; the Brans-Dickc 
parameter is uj > 40000. 

The proposed future missions aimed at improving 
the accuracy of 7 are: 

1. GP-B (geodetic precession) — 10~ 5 . 

2. Bepi-Colombo (retardation) — 10~ 6 . 

3. GAIA (deflection) - (10~ 5 - 10~ 7 ). 

4. ASTROD I (W.-T. Ni) (retardation) — 10" 7 . 

5. LATOR (Turyshev et al.) — 10~ 8 . 

6. ASTROD (W.-T. Ni) — 10" 9 . 

There are three problems related to G , whose origin 
is mainly related to unified models predictions: 

1) absolute G measurements, 

2) possible time variations of G , 

3) possible scale variations of G — non-Newtonian, or 

new interactions. 

Absolute measurements of G . There are many lab- 
oratory determinations of G with errors of the order 
10~ 3 , and only 4 have been on the level of 10 -4 in the 
80s. They are given in Table 1 (in 10~ n m 3 kg -1 s -2 ). 

Table 1 

1. Facy and Pontikis, France 1972 6,6714 ± 0.0006 

2. Sagitov et al., Russia 1979 6,6745 ± 0.0008 

3. Luther and Towler, USA 1982 6,6726 ± 0.0005 

4. Karagioz, Russia 1988 6,6731 ± 0.0004 

From this table it is evident that the first three ex- 
periments contradict each other (the results do not over- 
lap within their uncertainties). And only the fourth ex- 
periment is in accord with the third one. 

The official CODATA value of 1986 was 

G = (6, 67259 ± 0.00085) • 10~ n m 3 kg- 1 s~ 2 , (11) 
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based on Luther and Towler's determination. But af- 
ter precise measurements of G by different groups, the 
situation became more vague. 

As one may see from the Cavendish conference data 
of 1998 [55]) the results of 7 groups could agree with each 
other only on the level of 1CT 3 . So, CODATA adopted 
in 1999 

G = (6.673 ± 0.001) • 10~ n m 3 kg" 1 s~ 2 . (12) 

The most recent and precise G measurements do 
not agree, and some of them differ from the CODATA 
value of 1986. They are given in Table 2 (in 10" 11 
m 3 kg -1 s~ 2 ) 

But, from 2004 CODATA gives: 

G = 6.6742(10) • 10" n m 3 kg" 1 s~ 2 . (13) 

So we see that we are not too far (a little more than one 
order) from Cavendish who obtained a value of G at the 
level of 10 -2 two centuries ago. The situation with the 
measurement of the absolute value of G is really dif- 
ferent from that with atomic constants values and their 
uncertainties. This means that either the limit of ter- 
restrial accuracies of determining G has been reached 
or we have some new physics enterfering the measure- 
ment procedure [7 ■ The former means that maybe we 
should turn to space experiments to measure G [121111] , 
while the latter means that a more thorough study of 
theories generalizing Einstein's GR or unified theories is 
necessary. 

There also exist some satellite determinations of G 
(namely, G ■ AfEarth) at the level of 10 -9 (so, could we 
know G much better, our knowledge of masses of the 
Earth and other planets and consequently their models 
would be much better). There are also several less pre- 
cise geophysical determinations in mines at the level of 
10 -3 , but they do not improve the situation. 

A precise knowledge of G is necessary, above all, 
since it is a FPC; next, for evaluation of mass of the 
Earth, planets, their mean density and, finally, for con- 
struction of Earth models; for transition from mechani- 
cal to electromagnetic units and vice versa; for evalua- 
tion of other constants through relations between them 
obtained by unified theories; for finding new possible 
types of interactions and geophysical effects; for some 
practical applications like increasing of modern gra- 
diometer precision, since they demand calibration by a 
gravitational field of a standard body depending on G: 
the high accuracy of their calibration ( 10~ 5 - 10 -6 ) 
requires the same accuracy of G . 

3. The knowledge of constants values has not only a fun- 
damental but also a metrological meaning. The modern 
system of standards is mainly based on stable physical 
phenomena. So, the stability of constants plays a cru- 
cial role. And it may be even more important if new 
definitions of the units via fixed fundamental constants 
in the International System of Units (SI) will be intro- 
duced in 2011, as is suggested now. Let us give some 



short historical references to the evolution of systems of 
units. 

Before the official adoption of the Metric Conven- 
tion (1875) in 1799, two platinum basic standards (me- 
tre and kilogram) were introduced and placed in the 
Archive of the French Republic, which started the Met- 
ric System. 

Gauss in 1832 created the first coherent system of 
units introducing the second and measured magnetic 
fields in terms of mechanical units: mm, gram and sec- 
ond. 

In 1874, the British Association for Development of 
Science introduced the coherent system of units CGS 
based on Maxwell's and Thomson's suggestions of 1860: 
centimetre — gram — second. 

In 1880, the British Association and International 
Congress of Electricians adopted practical units — ohm, 
volt and ampere. 

In 1889, the First General Committee on Weights 
and Measures (BIPM) adopted the MKS System (metre, 
kilogram, second). 

In 1939, the MKSA system (MKS + ampere) was 
suggested (based on Georgi's suggestion of 1901) and 
officially adopted by BIPM in 1946. 

In 1954. kelvin and candela were introduced by 
BIPM, and the whole system of 6 basic units was called 
SI (International System of Units) in 1960. 

In 1971, the mole was added to SI, so up till now 
SI has 7 basic units: metre, kilogram, second, 
ampere, kelvin, candela and mole. 

A new stage of SI evolution started in October 
2005 when the International Committee for Weights 
and Measures (CIPM) adopted a recommendation on 
preparative steps towards redefining the 

kilogram, ampere, kelvin and mole, 
so that these units be linked to exactly known values of 
fundamental constants. Mills et al. 99J proposed that 
these four base units should be given new definitions 
linking them to exactly defined values of the Planck 
constant h, elementary charge e, the Boltzmann con- 
stant k and the Avogadro constant Na , respectively. 
This would mean that six of the seven base units of 
the SI would be defined in terms of true invariants of 
nature. Not only would these four fundamental con- 
stants have exactly defined values, but also the uncer- 
tainties of many of the other fundamental constants of 
physics would be either eliminated or appreciably re- 
duced. They suggested wordings for the four new def- 
initions and presented a novel way to define the entire 
SI units explicitly using such definitions without making 
any distinction between base units and derived units: 

The metre, unit of length, is such that the speed 
of light in vacuum c is exactly 

299 792 458 metres per second. 

Such a definition is simple, concise and makes clear 
the fundamental constant to which the unit is linked and 
the exact value of that constant. 
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Table 2 

1. Gundlach and Merkowitz (USA) [M] 

2. Armstrong (New Zealand, MSL) 

3. Karagioz (Moscow, Russia) 

4. Luo Zhun (Wuhan, China) 

5. Quinn ct al. (BIPM) (1) 

(2) 

6. Schlammingcr et al. (CH) 

If this general form were chosen, it would be appro- 
priate to choose definitions of the same form for all seven 
base units. 

Thus, for the second, candela, kilogram, ampere, 
kelvin and mole we would have: 

The second, unit of time, is such that the ground 
state hyperfine splitting transition frequency of the cae- 
sium 133 atom is exactly 
9 192 631 770 hertz. 

The candela, unit of luminous intensity in a 
given direction, is such that the spectral luminous effi- 
cacy of monochromatic radiation of frequency 540 x 1012 
hertz is exactly 

683 lumens per watt. 

The kilogram, unit of mass, is such that the Planck 
constant is exactly 

6.6260693 x 10~ 34 joule second. 

The ampere, unit of electric current, is such that 
the elementary charge is exactly 
1.60217653 x 10" 19 coulomb. 

The kelvin, unit of thermodynamic temperature, 
is such that the Boltzmann constant is exactly 

1.3806505 x 10~ 23 joule per kelvin. 

The mole, unit of amount of substance of a spec- 
ified elementary entity, which may be an atom, 
molecule, ion, electron, any other particle or a speci- 
fied group of such particles, is such that the Avogadro 
constant is exactly 

6.022 1 4 1 5 x 10 23 per mole. 

Of course, there still remain a lot of problems to be 
solved before real introduction of these new definitions 
(raising the precision of absolute values of some con- 
stants, choosing the variant of kg realization, testing 
the self-consistency of definitions, possible stability of 
constants, readiness of the world community to accept 
these changes (science, industry, trade, educational level 
etc.) 

4. Time Variations of G. As all physical laws were 
established and tested during the last 2-3 centuries in 
experiments on Earth and in the near space, i.e., at 
rather short space and time intervals as compared with 
the radius and age of the Universe, the possibility of 
slow variations of constants (i.e., with the rate of the 
evolution of the Universe or slower) cannot be excluded 



2000 6, 674215 ±0,000092 

2002 6, 6742 ±0,0007 

2003 6.6729 ±0.0005 
1998 6, 6699 ±0,0007 

2001 6, 6693 ± 0, 0009 

2001 6, 6689 ±0,0014 

2002 6.7404 ±0.000033 

a priori. So, the assumption of absolute stability of con- 
stants is an extrapolation, and each time we must test 
it. 

The problem of FPC variations arose with the at- 
tempts to explain the relations between micro- and 
macro-world phenomena. Dirac was the first to in- 
troduce (in 1937) the so-called "Large Number Hypoth- 
esis" which relates some known very big (or very small) 
numbers with the dimensionless age of the Universe 
T ~ 10 40 (age of the Universe in seconds 10 17 , divided 
by the characteristic elementary particle time 10 -23 sec- 
onds). He suggested that the ratio of the gravitational 
to strong interaction strengths, Gm^/(hc) ~ 10~ 40 , 
is inversely proportional to the age of the Universe: 
Gnip/hc ~ T _1 . Then, as the age varies, some con- 
stants or their combinations must vary as well. The 
atomic constants seemed to Dirac to be more stable, so 
he chose variation of G as T -1 . 

After the original Dirac hypothesis, some new ones 
appeared (Gamow, Teller, Landau, Terazawa, Staniuko- 
vich etc., see [4UT2] ). as well as some generalized theories 
of gravitation admitting variations of an effective gravi- 
tational coupling. We can single out three stages in the 
development of this field: 

1. Studies of theories and hypotheses with FPC vari- 
ations, their predictions and confrontation with 
experiments (1937-1977). 

2. Creation of theories admitting variations of an ef- 
fective gravitational constant in a particular sys- 
tem of units, analyses of experimental and obser- 
vational data within these theories [73113] (1977- 
present). 

3. Analyses of FPC variations within unified models 
0I5III] (present). 

In the development of the first stage from the anal- 
ysis of the whole set of the then available astronomi- 
cal, astrophysical, geophysical and laboratory data, a 
conclusion was made [73l[74] that variations of atomic 
constants were excluded at the level of 10 -15 per year, 
but variations of the effective gravitational constant in 
the atomic system of units do not contradict the ex- 
perimental data on the level of 10~ 12 -r- 10~ 13 year _1 . 
Moreover, in [73l[74j[75] , a conception was worked out 
that variations of constants are not absolute but de- 
pend on the system of measurements (choice of stan- 
dards, units and devices using this or that fundamen- 
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tal interaction). Each fundamental interaction, through 
dynamics described by the corresponding theory, defines 
the system of units and the corresponding system of ba- 
sic standards, e.g., atomic and gravitational (ephemeris) 
seconds. 

Earlier reviews of some hypotheses on variations of 
FPCs and experimental tests can be found in [H[M7]- 

5. There are different astronomical, geophysical and 
laboratory data on possible FPC variations [12] . 

Astrophysical data. Here follow some recent ones. 

Comparing the data from absorption lines of atomic 
and molecular transition spectra in high-redshift QSO's, 
Varshalovich and Potckhin (Russia) [55] obtained for 
z = 2.8-3.1: 

\a/a\ < 1.6- 1CT 14 year" 1 , (14) 
and Drinkwater et al. [83 : 

\a/a\ < 10" 15 year" 1 for z = 0.25 (15) 
and 

\a/a\ < 5 ■ 1CT 16 year" 1 for z = 0.68 (16) 

for a model with zero deceleration parameter and H = 
75km -s" 1 -Mpc" 1 . 

A less precise conclusion was made on the basis of 
geophysical data. Indeed, according to F. Dyson (1972), 
/3-decay of Rei^ to Osi87 give: 

\afa\ < 5 • 10" 15 year" 1 . (17) 

We must point out that all astronomical and geo- 
physical estimations are strongly model-dependent. So, 
of course, it is always desirable to have laboratory tests 
of FPC variations. 

a) Such a test was first done by a Russian group in 

the Committee for Standards (Kolosnitsyn, 1975). 
Comparing rates of two different types of clocks, 
one based on a Cs standard and another on a beam 
molecular generator, they found that 

|d/a| < 10" 10 year" 1 . (18) 

b) More recent data were obtained by J. Prestage et al. 

[84] by comparing mercury and i?-maser clocks. 
Their result is 

\a/a\ < 3.7- 10~ 14 year" 1 . (19) 

All these limits were placed on the fine structure con- 
stant variations. From the analysis of decay rates of K 40 
and Reisj , a limit on possible variations of the weak in- 
teraction constant was obtained (see ab approach for 
variations of (3, e.g., in [76 ) 

\$/p\ < 10" 10 year" 1 . (20) 



But the most stringent data on variations of strong 
(Gs), electromagnetic and week (Gw) interaction con- 
stants were obtained by A. Schlyachter (USSR) in 1976 
from an analysis of the ancient natural nuclear reactor 
data in Gabon, Oklo, because the event took place 2-10 9 
years ago. They are as follows: 

\Gs/G s \ <5-10" 19 year" 1 , 
\a/a\ < 10" 17 year" 1 , 
\Gp/Gp\ < 2 • 10" 12 year" 1 . (21) 

Some studies of the strong interaction constant and 
its dependance on transfered momenta may be found 
in [115j . A recent review on variations of a see in [97) . 

There have appeared some data on a possible varia- 
tion of a on the level of 10" 16 at some z [HE]- Other 
groups do not support these results. There also ap- 
peared data on possible variation of m e /m p (Var- 
shalovich et al.) The problem may be that even if they 
are correct, all these results are mean values of varia- 
tions at some epoch of the evolution of the Universe (a 
certain z interval). In essence, variations may be differ- 
ent at different epochs (if they exist at all), and at the 
next stage observational data should be analyzed with 
the account of evolution of the corresponding ("true"?) 
cosmological models. 

6. Now we are still having no unified theory of all four 
interactions. So it is possible to construct systems of 
measurements based on any of these four interactions. 
But in practice it is now done on the basis of the mostly 
worked out theory, electrodynamics (more precisely, on 
QED). Of course, it may also be done on the basis of the 
gravitational interaction (as was partly the case earlier). 
Then, different units of basic physical quantities arise, 
based on dynamics of a particular interaction, i.e., the 
atomic (electromagnetic) second, defined via frequency 
of atomic transitions, or the gravitational second de- 
fined by the mean Earth motion around the Sun (the 
ephemeris time). 

It does not follow from anything that these two sec- 
onds are always synchronized in time and space. So, 
in principle, they may evolve relative to each other, for 
example, at the rate of the evolution of the Universe or 
at some slower rate. 

That is why, in general, variations of the gravita- 
tional constant are possible in the atomic system of units 
(c, h, m are constant, Jordan frame) and masses of all 
particles — in the gravitational system of units ( G , H , c 
are constant by definition, Einstein frame). In practice, 
we can test only the first variant since the modern basic 
standards are defined in the atomic system of measure- 
ments. Possible variations of the FPCs must be tested 
experimentally, but for this purpose it is necessary to 
have the corresponding theories admitting such varia- 
tions and their certain effects. 

Mathematically, these systems of measurement may 
be realized as conformally related metric forms. Arbi- 
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trary conformal transformations give us a transition to 
an arbitrary system of measurements. 

We know that the scalar-tensor and multidimen- 
sional theories are the corresponding frameworks for 
these variations. So, one of the ways to describe vari- 
able gravitational coupling is the introduction of a scalar 
field as an additional variable of the gravitational inter- 
action. It may be done by different means (e.g. Jordan, 
Brans-Dicke, Canuto and others). We have suggested a 
variant of gravitational theory with a conformal scalar 
field (Higgs-type field [77ll4]). where Einstein's GR may 
be considered as a result of spontaneous breaking of 
the conformal symmetry (Domokos, 1976) [4]. In our 
variant, a spontaneous breaking of the global gauge in- 
variance leads to a nonsingular cosmology [75). Besides, 
we may get variations of the effective gravitational con- 
stant in the atomic system of units when m , c , h are 
constant and variations of all masses in the gravitational 
system of units (G, c, h are constant). It is done on 
the basis of approximate [79] and exact cosmological 
solutions with local inhomogenity [5U] , 

The effective gravitational constant is calculated us- 
ing the equations of motion. Post-Newtonian expansion 
is also used in order to confront the theory with ex- 
isting experimental data. Among the post-Newtonian 
parameters, the parameter / describing variations of G 
is included. It is defined as 



1 d{GM) 
GM Jt 



fH. 



(22) 



According to Hellings' data [81] from the Viking mission, 
7 - 1 = (-1.2 ± 1.6) • 10 3 , / = (4 ± 8) • 10~ 2 .(23) 



In the theory with a conformal Higgs field 79,80], we 
have obtained the following relation between / and 7 : 



/ = 4(7-1). 



(24) 



Using Hellings' data for 7, we can calculate / in our 
variant and compare it with / from [81] . Then we get 
/ = (-9, 6 ± 12, 8) • 10~ 3 , which agrees with (24) within 
its accuracy. 

We have used here only Hellings' data on variations 
of G . Other theoretical calculations in different models 
give the following predictions: less than 10~ 12 per year 
in multidimensional models 9,89,91,92], less than 10 -14 
per year [517] . 

But the situation with the experiment and observa- 
tions is not so simple. Along with [ST], there are some 
other data [4,5, 12 , but the most precise are: 



1. Hellings' result 

\G/G\ < (2 ±4) • 1CT 12 year" 1 . (25) 

2. A result from nucleosythesis (Acceta et al., 1992): 

\G/G\ < (±0.9) • 10~ 12 year" 1 . (26) 



3. E.V. Pitjeva's result (Russia) [S3], based on satel- 
lites and planets motion: 



\G/G\ < (0±2) • 10~ 12 year" 1 



(27) 



4. Some new results from pulsars and Big Bang nu- 
cleosynthesis (BBN) at the level of 10 -12 per year. 

There are also BBN data of Copi et al., 2003: 



3 • 10~ 13 < G/G < 4 - 10 



-13 



(28) 



As to other experimental or observational data, the 
results are of different quality. The most reliable ones 
are based on lunar laser ranging (LLR) (Muller et al, 
1993, Williams et al, 1996, Nordtvedt, 2003). They are 
not better than 10~ 12 per year. Here, once more we 
see that there is a need for corresponding theoretical 
and experinmental studies. Probably, future space mis- 
sions like Earth SEE-satellite fl ^ fTT | [T2] , [T3] or missions 
to other planets and lunar laser ranging will be a deci- 
sive step in solving the problem of temporal variations of 
G and determining the fates of different theories which 
predict them, since the greater is the time interval be- 
tween successive masurements and, of course, the more 
precise they are, the more stringent results will be ob- 
tained. 

As we saw, different theoretical schemes lead to tem- 
poral variations of the effective gravitational constant: 

1 . Empirical models and theories of Dirac type, where 
G is replaced with G{t) . 

2. Numerous scalar-tensor theories of Jordan-Brans- 
Dicke type where G depends on the scalar field 
a(t). 

3. Gravitational theories with a conformal scalar field 
arising in different approaches [73l[74l, Hl [77[[56] . 

4. Multidimensional unified theories in which there 
are dilaton fields and effective scalar fields appear- 
ing in our 4-dimensional spacetime from additional 
dimensions [87L lTj . They also may help in solving 
the problem of a variable cosmological constant 
changing from very large Planckian to present (0.7 
of total energy) values. 

As was shown in [S][S7l[T], temporal variations of 
FPCs are connected with each other in multidimensional 
models of unification of interactions. So, experimental 
tests on a/ a may at the same time be used for esti- 
mation of G/G and vice versa. Moreover, variations of 
G are also related to the cosmological parameters p , SI 
and q , which gives opportunities of raising the precision 
of their determination. 

Since FPC variations are closely connected with the 
behaviour of internal scale factors, it is a direct probe 
of properties of extra dimensions and the corresponding 
unified theories [H11H1H] • From this point of view, it is an 
additional test of not only gravity and cosmology, but 
of unified theories of physical interactions as well. 
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7. Non-Newtonian interactions, or range varia- 
tions of G . Nearly all modified theories of gravity and 
unified theories also predict some deviations from the 
Newton law (inverse square law, ISL) or composition- 
dependent violations of the Equivalence Principle (EP) 
due to appearance of new possible massive particles 
(partners) [5 . Experimental data exclude the existence 
of these particles on a very good level in nearly all ranges 
except less than the millimetre and metres and hundreds 
of metres ranges. Our recent analysis of experimental 
bounds and new limits on possible ISL violation using 
a new method and modern precession data from satel- 
lites, planets, binary pulsar and LLR data were obtained 
in [!55]. 

In the Einstein theory G is a true constant. But, 
if we think that G may vary with time, then, from a 
relativistic point of view, it may vary with distance as 
well. In GR massless gravitons are mediators of the 
gravitational interaction, they obey second-order differ- 
ential equations and interact with matter with a con- 
stant strength G. If any of these requirements is vio- 
lated, we come in general to deviations from the Newton 
law with range (or to a generalization of GR). 

In [5] we have analyzed several classes of such theo- 
ries: 

1. Theories with massive gravitons like bimetric ones 
or theories with a A-term. 

2. Theories with an effective gravitational constant 
like the general scalar-tensor ones. 

3. Theories with torsion. 

4. Theories with higher derivatives (4th-order equa- 
tions etc.), where massive modes appear, leading 
to short-range additional forces. 

5. More elaborated theories with other mediators in 
addition to gravitons (partners), like supergravity, 
superstrings, M-theory etc. 

6. Theories with nonlinearities induced by any known 
physical interactions (Born-Infeld etc.) 

7. Phenomenological models where a detailed mech- 
anism of deviation is not known (fifth or other 
force) . 

8. Modifications of the Newton law at large ranges 
(MOND etc.), small acceleration at a > a,Q (Pio- 
neer anomaly, etc.) 

In all these theories, some effective or real masses appear 
leading to Yukawa-type (or power-law) deviations from 
the Newton law, characterized by the strength a and 
the range A . 

There exist some mo del- dependent estimates of these 
forces. The most well-known one belongs to Scherk 
(1979) from supergravity where the graviton is accom- 
panied by a spin-1 partner (graviphoton) leading to an 
additional repulsion. Other models were suggested by 
Moody and Wilczek (1984) — introduction of a pseudo- 
scalar particle leading to an additional attraction be- 
tween macro-bodies with the range 2-10 -4 cm < A < 20 



cm and strength a from 1 to 10~ 10 in this range. An- 
other supersymmetric model was elaborated by Fayet 
(1986, 1990), where a spin-1 partner of a massive gravi- 
ton gives an additional repulsion in the range of the 
order 10 3 km and a of the order 10 -13 . 

A scalar field to adjust A was introduced by S. Wein- 
berg in 1989, with a mass smaller than 10 _3 eV/c 2 , or a 
range greater than 0.1 mm. One more variant was sug- 
gested by Peccei, Sola and Wetterich (1987), leading to 
additional attraction with a range smaller than 10 km. 

Some p-brane models (ADD, brane worlds) also 
predict non-Newtonian additional interactions of both 
Yukawa or power-law types, in particular, in the sub- 
millimetre range, which is intensively discussed nowa- 
days 13,96]. On PPN parameters for multidimensional 
models with p-branes see above, Sec. 2. 

The Pioneer anomaly. A more serious evidence on a 
possible violation of Newton's law has come to us from 
space, namely, from data processing on the motion of the 
spacecrafts Pioneer 10 and 11, referring to length ranges 
of the order of or exceeding the size of the Solar system. 
The discovered anomalous (additional) acceleration is 

M\ 

(8.60 ± 1.34)-10- 8 cm/s 2 , 

it acts on the spacecrafts and is directed towards the 
Sun. This acceleration is not explained by any known 
effects, bodies or influences related to the design of the 
spacecrafts themselves (leakage etc.), as was confirmed 
by independent calculations. 

Many different approaches have been analyzed both 
in the framework of standard theories and invoking 
new physics (Chongming Xu's talk at ICGA-7, Taiwan, 
November 2005), but none of them now seems to be 
sufficiently convincing and generally accepted. There 
are the following approaches using standard physics: 

— an unknown mass distribution in the Solar system 

(Kuiper's belt), interplanetary or interstellar dust, 
local effects due to the Universe expansion [100 ; 

— employing the Schwarzschild solution with an ex- 

panding boundary |101U102j etc. 

Among the approaches using new physics one can 
mention: 

— a variable cosmological constant [103] ; 

— a variable gravitational constant |106j ; 

— a new PPN theory connecting local scales with the 

cosmological expansion |104j ; 

— the five-dimensional Kaluza-Klein (KK) theory with 

a time- variable fifth dimension and varying funda- 
mental physical constants [105] : 

— Moffat's |107j non-symmetric gravitational theory; 

— Milgrom's [108, 109 modified Newtonian dynamics 

(MOND); 

— special scalar-tensor theories of gravity [110] ; 
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— approaches using some ideas of multidimensional 

theories [111] ; 

— modified general relativity with a generalized stress- 

energy tensor [1 12] etc. 

This Pioneer anomaly has caused new proposals of 
space missions with more precise experiments and a wide 
spectrum of research at the Solar system length range 
and beyond: 

— Cosmic Vision 2015-2025, suggested by the Euro- 

pean Space Agency, and 

— Pioneer Anomaly Explorer, suggested by NASA 

urn 

So, we hope they can contribute a lot to our knowl- 
edge of gravity and unified models. 

8. Space project SEE and laboratory projects. 

We have seen that there are three problems connected 
with G . There is a promising new multi-purpose space 
experiment SEE (Satellite Energy Exchange) flQlITT] 
which addresses all these problems and may be more 
effective in solving them than other laboratory or space 
experiments. 

This experiment is based on a limited 3-body prob- 
lem of celestial mechanics: small and large masses in a 
drag-free satellite and the Earth. Unique horse-shoe or- 
bits, which are effectively one-dimensional, are used in 
it. 

The aims of the SEE-project are to measure the In- 
verse Square Law (ISL) and the Equivalence Principle 
(EP) at ranges of metres and the Earth radius, G and 
the absolute value of G with unprecedented accuracies. 

We have studied many aspects of the SEE project 
[THE]: 

1 . A wide range of trajectories with the aim of finding 
optimal ones: 

— circular in a spherical field; 

— the same plus the Earth's quadrupole modes; 

— elliptic with eccentricity less than 0.05. 

2. Estimations of other celestial bodies' influence. 

3. Estimation of the sensitivity of trajectories to 
changes in G and the Yukawa-type interaction 
strength parameter a . 

4. Modelling the measurement procedures for G and 
a by different methods for different ranges and for 
different satellite altitudes: 

— optimal, 1500 km, 

— ISS free flying platform, 500 km and 

— a high orbit, 3000 km. 

5. Estimations of some sources of error were made 
for: 

— radial oscillations of the shepherd's surface; 

— longitudeal oscillations of the capsule; 



— transversal oscillations of the capsule; 

— shepherd's non-sphericity; 

— limits on the quadrupole moment of the shep- 

herd; 

— limits on admissible charges and time scales of 

charging by high energy particles etc. 

6. Error budgets for G, G and G(r) were calculated. 

On the basis of all these studies, the general conclusion 
was that realization of the SEE project may improve our 
knowledge of G, G and G{r) by 3-4 orders of magni- 
tude. 

A laboratory experiment was also suggested in our 
paper [114] to test the possible range variations of G . 
It is an experiment for possible detection of new forces, 
or test of the inverse square law, parameterized by a 
Yukawa- type potential with strength a and range A. 
The installation comprises a ball with a spherical cav- 
ity whose centre is shifted with respect to the ball cen- 
tre. The ball is placed on a turn-table being subject to 
uniform rotation. A torsion balance as a sensitive ele- 
ment is placed inside the cavity. A uniform gravitational 
field created inside the ball does not affect the balance, 
but any non-gravitational forces create a torque which 
acts periodically during the rotation of the ball. The 
spectrum of harmonics was calculated. It is shown that 
preferable to use is the first harmonic in the measure- 
ments. Sensitivity of the method was evaluated, which 
is limited by uncertainties due to manufacturing of the 
elements and temperature fluctuations of the sensitive 
element. It was shown that the sensitivity of the method 
suggested may be at the level of 10 -10 in a in the range 
of A from 0.1 to 10 7 m in the space of the Yukawa pa- 
rameters (a, A) . 
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